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The  exceedance  point  process  approach  of 
Hsing  et  al.  is  extended  to  multivariate 
stationary  sequences  and  some  weak  con- 
vergence results  are  obtained.  It  is  well 
known  that  under  general  mixing  assump- 
tions, high  level  exceedances  typically 
have  a limiting  Compound  Poisson  structure 
where  multiple  events  are  caused  by  the 
clustering  of  exceedances.  In  this  paper  we 
explore  (a)  the  precise  effect  of  such 
clustering  on  the  limit,  and  (b)  the  relation- 
ship between  point  process  convergence 
and  the  limiting  behavior  of  maxima.  Fol- 


lowing this,  the  notion  of  multivariate  ex- 
tremal index  is  introduced  which  is  shown 
to  have  properties  analogous  to  its  uni- 
variate counterpart.  Two  examples  of 
bivariate  moving  average  sequences  are 
presented  for  which  the  extremal  index  is 
calculated  in  some  special  cases. 
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1.  Introduction 

Extreme  value  theory  for  multivariate  iid  sequences 
has  been  studied  for  quite  some  time  now  but  attention 
to  the  dependent  case  has  been  relatively  recent.  For 
univariate  sequences  it  is  known  that  local  dependence 
causes  extreme  values  to  occur  in  clusters,  which  in  turn 
results  in  a stochastically  smaller  distribution  for  the 
maximum  than  if  the  observations  were  independent.  We 
begin  with  a brief  review  of  these  results,  which  we 
shall  later  extend  to  the  multivariate  case. 

Let  {£„ } be  a univariate  stationary  sequence.  Write 
M„=max{^l,...,^n } and  for  t>0,  let  {«„( r)}  denote  a 
sequence  satisfying  linv^nP  {£,>«„  (t)}=t.  Under 
quite  general  mixing  assumptions  there  exist  constants 
O<0'<0"<1  such  that 

lim  sup  P {M„^Mn(r)}=e”fl’T  and 

n— 

lim  inf  P {Mn<un(r)}=eirT 


for  all  r.  (See  Ref.  [1],  although  the  idea  actually  dates 
back  to  Refs.  [2-4].)  Thus  if  P {Mn^un(T0)}  converges 
for  some  t0,  then  8'=8'\=8,  say)  and  hence 
lim M®/){iW„s«„(T)}=e‘'r  for  all  r>0.  The  common 
value  8 is  then  called  the  extremal  index  of  {£„  }•  We 
shall  assume  8 to  be  positive  whenever  it  exists,  since  the 
case  8=0  corresponds  to  a degenerate  limiting  distribu- 
tion for  Mn . Note  that  8=  1 for  iid  sequences.  Let  {£„  } be 
an  iid  sequence  with  £i=d  called  the  associated  iid 
sequence , and  write  M„= max{£ }.  If  {£„  } has  ex- 
tremal index  8 and  lim,,.,,*  P {Mn<vn(t)}=H (t)  for  a 
suitable  family  of  normalizing  constants  {v„(0},  then  it 
follows  (upon  identifying  e ~6t  with  //(/))  that  limn^ 
P {Mn<vn(t)}=H{t)  where 

H(t)=H(t)e.  (1) 

The  extremal  index  is  thus  a measure  of  the  effect  of 


543 


Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


dependence  on  the  limiting  distribution  of  Mn.  The 
stochastically  smaller  limiting  distribution  of  M„  is  in 
fact  a direct  result  of  the  clustering  of  extremes,  as 
explained  below.  See  Ref.  [5]  for  details. 

For  fixed  r> 0 let  the  exceedance  point  process 
Nn=N(nT)  be  defined  by 

X„{B  )=X  I {Hn&Bh  5C  [0,1], 

/=■  1 

where  lA  denotes  the  indicator  function  of  the  event  A . 
Then  for  a broad  class  of  weakly  dependent  sequences, 
the  limit  in  distribution  of  N„,  if  it  exists,  is  a Compound 
Poisson  process  with  intensity  0rand  multiplicity  distri- 
bution TTon  {1,2,...}.  The  Poisson  events  may  in  fact  be 
regarded  as  the  positions  of  “exceedance  clusters” 
while  the  multiplicities  correspond  to  cluster  sizes. 
More  explicitly,  one  may  divide  the  n observations  into 
kn  blocks  of  roughly  equal  size  and  regard  exceedances 
within  each  block  as  forming  a single  “cluster”,  so  that 
the  cluster  sizes  are  given  by  Nn(Jt),  i=\,...,kn,  where 
J X=J n.i)=(n~  , *-]•  For  a suitable  choice  of  kn  depending 
on  the  mixing  rate  of  {£„ },  one  then  has 

lim  P {Nn  (/,)=/  | Nn(J{)>0}=Tr(j),  1, 

n — >°° 

and 

lim  P‘-{y,(Jl)=0}=liml,R[0,l]=0} 

/J—>CO  AI— >00 

=lim  P {Mn^un(T)y=e~eT, 

/?— >oo 

so  that  in  particular,  lim„^  kn  P {Nn(J\)>0}=07.  Hence 
lim  E A„[0,l]=lim  knENn{J ,) 

n—>  co  n — >oo 

=lim  knE (N„ (./ , ) | Nn(Ji>0)P  {Nn(J {)>()} 

n— >oo 

=6t  lim  E(Nn(J[)\Nn(Jl)>0), 

n— >oo 

while  ori  the  other  hand,  lim^a,  L/VTCU^lim,,,.^ 
nP  {^i>«„(r)}=r.  The  cluster  size  distribution  and  the 
extremal  index  are  therefore  related  by 

lim  E(Nn(Jf  | A„(7i)>O)=l/0.  (2) 

n— >oo 

Now  let  nE:2Z}  be  a multivariate  sta- 

tionary sequence  where  zfe  1 is  a fixed  integer,  and 
write  where  M„rmax{£,, }, 

j=\,...,d.  The  study  of  multivariate  extremes  began  in 
the  early  1950s,  focusing  mainly  on  the  limiting  behav- 
ior of  M„  under  a linear  normalization,  when  the  obser- 
vations are  iid.  The  resulting  class  of  limiting  distribu- 
tions was  characterized  in  Ref.  [6]  and  domains  of 
attraction  criteria  were  given  in  Ref.  [7],  See  also  Ref. 


[8],  Chapter  5,  for  an  account  of  the  literature  surround- 
ing this  theory.  For  stationary  sequences  satisfying  a 
general  mixing  assumption,  it  is  known  (see  Refs.  [9, 
10],  and  Theorem  1.1  below)  that  the  class  of  limiting 
distributions  of  Mn  is  the  same  as  for  iid  sequences.  In 
this  paper  we  explore  the  precise  effect  of  dependence 
on  the  limiting  distribution  by  extending  the  univariate 
theory  described  above  to  the  multivariate  case.  Essen- 
tially, this  involves  studying  the  inter-relationship  be- 
tween the  two  dependence  structures  present,  one  due  to 
dependence  over  time  and  the  other  due  to  the  depen- 
dence between  the  various  components  of  the  multivari- 
ate observations.  The  ideas  become  most  transparent 
when  presented  in  terms  of  so-called  dependence  func- 
tions [8],  Here  we  adopt  the  slightly  modified  definition 
found  in  Ref.  [9].  A distribution  function  D on  [0,1  ]rf  is 
called  a dependence  function  if  Dj(Dj(u))=Dj(u), 
wE[0,l],y=l,...,<7,  where  the  subscript  j signifies  the  j th 
marginal.  The  dependence  function  of  a distribution  F 
on  IRd  is  defined  by  DF(u)=P  {FfXA^ 
uu.-.J7d(Xd)^ud},  u=(ui,...,ud)E  [0,1  ]d,  where 

(X|,...,Xrf)  is  a random  vector  with  distribution  F . More 
generally,  any  dependence  function  satisfying 
F (x)=D  (Fi(x[),...,Fd(xd))  could  be  defined  to  be  a de- 
pendence function  of  F , although  the  present  choice  is 
a natural  one. 

Write  r=(0,l)A {1  > where  l=(l,...,l)ElRrf,  and  for 
t=(r,,...,r(/)E7\  let  vn(t)=(v„i(ti),...,vnd(td))  where  vnj(tj) 
satisfies  lim„_^  nP  { p>vnj (tj)}=~ logf; . Let  H„  denote 
the  distribution  function  of  Mn  (i.e.,  H n(x)=P  {M n^x  }), 
with  marginals  Hnj,  j=\,...,d.  Then  (see  Refs.  [8,  11]), 

//„(v„(t))->"//(t)  (3) 

if  and  only  if 

kf„j(vnj(tj))-^>"  ffftj),  j=\ ,...,d , and  DHn{ t)— >“  DH{ t). 

The  limiting  behavior  of  Mn  can  therefore  be  separated 
into  two  parts,  one  pertaining  to  the  convergence  of  the 
marginals  (a  univariate  problem)  and  the  other  to  the 
convergence  of  the  dependence  functions.  Here  we  fo- 
cus attention  exclusively  on  the  latter.  It  should  be  noted 
that  the  choice  of  normalising  constants  does  not  affect 
the  dependence  function  of  the  limit  distribution  H , but 
only  alters  the  marginals  (see  Ref.  [9],  Lemma  3.2). 
Since  our  main  interest  is  in  the  dependence  function, 
the  present  choice  of  normalising  constants  is  appropri- 
ate in  view  of  the  fact  that  it  results  in  Uniform[0,l] 
marginals  for  the  limit  distribution  when  {£„}  is  iid,  so 
that  in  particular  Dh=H . According  to  Theorem  3.3  of 
Ref.  [9],  the  class  of  all  possible  limits  H in  Eq.  (3)  (for 
iid  {f, })  is  precisely  the  class  of  extreme  dependence 
functions,  that  is  those  that  satisfy 
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Dn(t)=D(t",...,tj)  (4) 

for  each  «>1  and  t=(fi,...,/rf)E[0,l]d.  Theorem  1.1  below 
shows  that  the  same  is  true  also  if  {£„ } is  a stationary 
sequence  satisfying  the  following  mixing  condition. 

For  ter,  let 

^*(v«(t))=cr{(^>vn;(t;)) : k<i<l,  j=\,...,d} 
and  for  1 </<«—!,  define 

a„.,=sup{  \P(ADB)-P{A)P{B)  \ : AE^(vn( t)), 

BE§;+/(v„( t)),  1 <k<k+l<n  }. 

The  mixing  condition  A(v„(t))  is  then  said  to  hold  if 
a„jn— >0  for  some  sequence  {/„ } satisfying  /„//?— > 0.  This 
is  the  multivariate  version  of  the  mixing  condition  used 
in  Ref.  [5]  and  is  slightly  stronger  than  the  D(un)  condi- 
tion in  Ref.  [9].  Henceforth  {£„}  will  be  assumed  to 
satisfy  A(v„(t)),  for  some  or  all  t,  as  required. 

Theorem  1.1.  Let  {£„  } satisfy  A(v„(t))  for  all  t ET and 
suppose  that  P{M„<v„( t)}— H( t),  non-degenerate. 
Then  DH  is  an  extreme  dependence  function  and  hence, 
in  particular,  H(tc)=Hc(t)  for  each  te[0,l]rf  and  c>0 
(where  t 

Proof:  The  first  part  is  an  immediate  consequence  of 
Theorem  4.2  of  Ref.  [9]  while  the  second  part  follows 
from  the  definition  of  extreme  dependence  functions 
upon  noting  that  (by  the  univariate  theory  described 
above),  the  marginals  of  H are  of  the  form  Hj(tj)=tf 
where  0,  is  the  extremal  index  of  { },  the  yth-compo- 
nent  sequence  of  {£„  }• 

In  the  next  section  we  apply  the  exceedance  point 
process  approach  to  multivariate  extremes  and  obtain 
some  weak  convergence  results.  The  multivariate  ex- 
tremal index  is  then  defined  (in  Sec.  3),  based  on  the 
multivariate  analogue  of  Eq.  (1).  It  is  seen  to  be  a func- 
tion of  only  d—  1 variables  and  its  properties  naturally 
extend  those  of  the  univariate  extremal  index.  Finally  in 
Sec.  4 we  consider  two  examples  of  bivariate  moving 
average  sequences  for  which  the  computation  of  the 
extremal  index  is  demonstrated. 

2.  Exceedance  Point  Processes 

Fix  tET and  let  6//=/{^>»n/^)>,  i=\,...,n,j=\,...d,  and  put 
8j=(8n,...,8id).  The  multivariate  exceedance  point  pro- 
cess Nn=N{nl)  is  then  defined  by 


M,(A)=S  BE[0,\].  (5) 

/-I 

Assume  that  {£„ } satisfies  A(v„(t)).  If  also  /V„— No, 
then  it  may  be  shown  (as  in  the  univariate  case)  that  the 
limit  N0  is  a point  process  on  [0,1]  which  is  of  Com- 
pound Poisson  type.  More  precisely,  the  Laplace  Trans- 
form of  N0  is  given  by 

-log£  exp{~X  f /,  dN0j}=v  f [ 

j~  1 *M0,1]  J xE\0J  r yGz^' 

d 

(l-exp{-X  y/(.v)»d7r(y)<iv.  (6) 

>1 

Here  N0j  denotes  the  yth-component  of  N0,  v is  a 
positive  constant,  77  is  a probability  distribution  on 
Z+'={0,1,2,...}‘\{0}  and^’s  are  non-negative  functions 
on  [0,1]. 

Let  {kn  } be  any  sequence  of  positive  integers  satisfy- 
ing 

kn — k„ln/n — >0,  and  knocnjn — ^0,  as  n — (7) 

Set  r„=[n/k„ ] (the  largest  integer  not  exceeding  n/kn) 
and  put  y„:=[0,r„/«  ].  Define  the  probability  distribution 
7 r„  on  TZf  by 

Trn(y)=P{Nn(Jnt)=y  I MTTMAO},  yETZf. 

The  following  theorem  which  gives  a useful  characteri- 
zation of  the  convergence  of  Nn  is  an  immediate  conse- 
quence of  the  results  in  Sec.  5 of  Ref.  [12]. 

Theorem  2.1.  N0  if  and  only  if  77„-ah'  77  and 

P {M„<  v„  (t) }— > e ",  and  in  that  case  the  Laplace  Trans- 
form of  No  is  given  by  Eq.  (6). 

Next  we  consider  the  iid  case  in  some  detail  and 
obtain  an  interesting  connection  with  Theorem  5.3.1  of 
Ref.  [8], 

Proposition  2.2.  Let  {£„}  be  iid  and  for  fixed  tET  let 
N„  be  defined  by  Eq.  (5).  If  N„—^d  No  then  the  multiplic- 
ity distribution  77  in  Eq.  (6)  is  supported  on  the  set 
S={0,1}'\{0>. 

Proof:  Observe  that  A(v„(t))  is  trivially  satisfied  since 
a„i=0  so  that  we  may  take  /„=  1 and  k„=n.  Then 
ir„(y)=P  {8\=y  I ^,^vn(t)},  yEZZf,  which  is  clearly  sup- 
ported on  S.  The  result  is  now  immediate  since  S'  is  a 
closed  set  and  77n— 77  by  Theorem  2.1.  □ 

Making  the  dependence  on  t explicit,  we  now  write 
Nn=N^ , N0=Nf\  v=vw  and  tt=ttw.  In  addition  we  shall 
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require  the  following  notation  from  Ref.  [8],  For 
1 <k<d.  let  j (k)=(jx,...Jk)  denote  a vector  with  integer- 
valued components  \<jx<ji<—<jk—d,  and  for 
x=(xx,...jcd)E.\RJ  write  Xj (*)=(*;>•••**/*)•  Define  the  “sur- 
vival function” 

G(x)=P  {£,  i>xu...Sd>xd  } 

and  write  Gm(x)=P  {^\jx>Xj &jk>xJk }.  For  each  j(&), 
let  yj(t)  denote  the  element  in  S={0,1  }A{0}  whose  j th 

component  equals  1 if  and  only  if  j=j,  for  some  j=l k. 

(This  defines  a natural  1-1  correspondence  between  S 
and  the  j(&)’s.) 

Theorem  2.3.  Let  {£„ } be  iid.  Then  Nf  for  some 
fixed  tG  T if  and  only  if 

lim  nGiik)(vn(t))=hm(t)<cc  (8) 

n— 

for  each  i(k),\—k^d.  In  that  case  has  Laplace 
Transform  given  by  Eq.  (6)  with 

d 

^(t)=Z(-l)*+l  2 *J(*>(‘)  (9) 

i-1  \sj<.<jksd 

and  with  7r(,)  determined  by  the  relations 

hm( t)=vw  2 ^(,>(y)-  (10) 

y-yj«  i 

Proof:  Write  Sk(t)=Hxsjt<...<j^dGm(t),  so  that 

vn (t) }=21i( ~ 1 )k+iSk (v„ (t)).  If  Eq.  (8)  holds  for  each  k, 

then 

d 

limnP{^^v„(t)}=2(-l)*+1  S (11) 

" k-i  1 ==/,<-</*<</ 

and  hence  limn^^/>{M„^v„(t)}=e_,<". 

Next  observe  that  for  each  j(&),  1 <k<d, 

Giik)(vn(t))=  2 P{8i=y}=P{^vn(t)}  2 7r„(y),(12) 

y-yj(o 

where  7ra(y)=P  {Sx=y  | ^i^v„(t)}.  Moreover,  this  rela- 
tionship is  invertible  in  the  sense  that  each  of  the  prob- 
abilities 77„(y),  y GS,  can  be  expressed  as  a linear  com- 
bination of  the  Gj(t)(v„(t))’s.  Therefore  by  Eq.  (8), 
lim„^7r„(y)=7r(,,(y)  (say)  exists  and  satisfies  Eq.  (10). 
Hence  by  Theorem  2.1  A^l)— A,1,11  where  /V,1,'1  has  the 
specified  parameters.  Conversely  if  N™— (Vo"  then  i t„ 
and  P{A/„<v„( t)}  converge  (by  Theorem  2.1  again), 
and  hence  Eq.  (8)  follows  by  virtue  of  Eqs.  (11)  and 
(12).  □ 


Corollary  2.4.  Let  {£„}  be  iid.  Then  Nf  -V  Nf  for 
each  t GT  if  and  only  if  P {Mn<vn{ t)}— >w  H{ t).  More- 
over H and  {vn\TT{"}teT  determine  each  other. 

Proof:  (Sketch)  The  first  part  follows  from  Theorem 
2.3  above  and  Theorem  5.3.1  of  Ref.  [8]  which  states 
that  P {M„<v„(t)}— >H’  H{ t)  if  and  only  if  Eq.  (8)  holds 
for  each  tG7\  Note  that  so  that  H and  the 

f("’s  can  be  obtained  from  each  other.  Also  the 
can  be  obtained  from  the  vw's  by  first  inverting  Eq.  (9) 
to  get  the  /2j(t,(t)’s  and  then  inverting  Eq.  (10).  (The 
inversion  of  Eq.  (9)  is  carried  out  inductively  using  the 
fact  that  the  weak  convergence  of  Hn(vn( t))  implies  that 
of  all  lower  dimensional  marginals.)  □ 

Analogous  results  for  the  dependent  case  take  on  a 
slightly  different  form.  Let  {£„ } be  a stationary  se- 
quence satisfying  A(v„(t))  for  each  tG7\  As  before  let 
rn=[n/kn]  where  {kn  } is  any  sequence  satisfying  Eq.  (7), 
and  define 

Gr,j(k)(vn  (t ))=P  x>V„j {(tj l),...JMrnjl>Vnjk(tjk  ) }. 

Theorem  2.5.  Let  {£„  } be  a stationary  sequence  satis- 
fying A(v„(t ))for  each  t E.T.  Then  P {MnSvn(t)}— H(t) 
if  and  only  if 

lim  knGrni(k](vn(t))=hm(t)<x 

n—>°° 

for  each  j(&),  1 and  tE7^  and  in  that  case 

d 

H{  t)=exp{2(-D*  S 

k-i  lsjt<—<jk£d 

Proof:  Observe  that  the  mixing  condition  A(v„(t))  im- 
plies that  {£„j(*)}  satisfies  A(v„j(i)(t))  for  each  j(fc)  (with 
obvious  notation).  Hence  it  may  be  shown  as  in  the 
univariate  case  (see  Lemma  2.1  of  Ref.  [1])  that 

P {Mnj{k)— vn^k)(i)}—Pk"  {A/r,j(*)— v„j(*)(t)}— : >0,  (13) 

for  each  j(A).  The  result  may  therefore  be  proved  in 
exactly  the  same  way  as  Theorem  5.3.1  of  Ref.  [8],  □ 

Remark:  Under  the  hypothesis  of  Theorem  2.5,  if 
Nf—yd  Nf\  t G7\  with  parameters  vw  and  7r(t)  then 
P{A7n<v„(t)}— >H  H( t)=e  as  in  the  iid  case.  However 
it  is  not  possible  in  general  to  recover  the  ^“’’s  from  H 
since  the  clustering  of  exceedances  may  cause  the  sup- 
port of  7T(t)  to  extend  beyond  S.  References  [9,  10]  give 
sufficient  conditions  (analogous  to  the  D\un)  condition 
of  Ref.  [13])  under  which  clustering  does  not  occur,  so 
that  Corollary  2.4  can  be  extended  to  stationary  se- 
quences satisfying  this  condition. 


546 


Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


A distribution  function  F on  IR''  is  said  to  be  indepen- 
dent if  F(xHI£,F;(x;),  xEIR^.  If  {£„}  is  iid  and 
P {A/„:£v„(t)}— >"  H( t),  then  it  follows  from  Corollary 
5.3.1  of  Ref.  [8]  that  H is  independent  if  and  only  if  the 
marginals  of  H are  pairwise  independent.  The  analogous 
result  for  the  dependent  case  is  stated  below.  The  proof 
(which  is  omitted)  is  essentially  the  same  as  for  the  iid 
case,  but  uses  Theorem  2.5  instead  of  Theorem  5.3.1  of 
Ref.  [8], 

Corollary  2.6.  Let  {£„  } be  a stationary  sequence  sat- 
isfying A(v„(t))  for  each  t E.T  and  suppose  that 
P{Mn^vn( t)}— >K’  H( t).  Then  H is  independent  if  and 
only  if  knP {Mrnj>vnj(tj),  Mrj>vnJ(t,)}-> 0 for  each 
l—j<l—d,  te  T,  i.e.,  if  and  only  if  knGr^k)(vn  (t)> — >0  for 
each  j(2)  and  each  t£7? 

It  is  shown  in  Ref.  [14]  that  H is  independent  if 
H{i)=Wdj.\Hj{tj)  for  some  tE(0,l)d.  Although  the  result 
in  [14]  only  stated  for  iid  sequences  under  a linear  nor- 
malization, the  proof  essentially  rests  on  the  defining 
property  of  extreme  dependence  functions,  namely  Eq. 
(1).  Consequently  the  result  extends  to  the  present  more 
general  situation  allowing  dependence  and  non-linear 
normalizations.  Corollary  2.6  can  therefore  be  improved 
as  follows. 

Corollary  2.7.  Let  {£„ } be  as  in  Corollary  2.6  and 
suppose  that  P{M„<v„(t)}— H (t).  Then  the  following 
are  equivalent: 

(i)  H is  independent, 

(ii)  H(X)=YidjmXHj(tj)  for  some  tE^l)'', 

(Hi)  knGr„nk)(vn (t)) — >0  for  each  j(2 ),for  some  tG(0,l)d. 

It  should  be  noted  that  Refs.  [9,  10]  give  some  inter- 
esting sufficient  conditions  for  H to  be  independent 
when  {£„ } is  a stationary  sequence.  A natural  question 
to  ask  in  the  present  context  is  whether  H is  independent 
whenever  H is.  Proposition  3.4  gives  a necessary  and 
sufficient  condition  for  this  in  terms  of  the  extremal 
index,  but  the  answer  in  general  is  negative  and  a coun- 
ter-example can  be  found  in  [10].  It  seems  more  plausi- 
ble that  the  converse  may  be  true,  i.e.,  that  H is  indepen- 
dent whenever  H is.  In  fact  however,  this  too  is  not  the 
case,  as  shown  by  an  interesting  counter-example  in 
[15]. 

We  conclude  this  section  by  stating  a result  which 
extends  Theorem  5.1  of  [5]  and  is  proved  similarly. 

Theorem  2.8.  Let  {f, } be  a stationary  sequence  satis- 
fying A(v„(t ))for  each  t£T and  suppose  that  Nf  — Nq] 
for  some  t GT  Then  Nf 1 for  each  c>0  and  fur- 

thermore, v{tC)=cvw  and  7r(‘C)=7r(,)  (where  tc=(t(,...,tj)). 


3.  The  Multivariate  Extremal  Index 

Let  {£„  } be  a stationary  sequence  and  {£„  } the  associ- 
ated iid  sequence.  Suppose  that  P {M„<v„(t )}—>“’  H( t) 
and  P {Mn<v„ (t )}—>"'  H( t).  The  multivariate  extremal 
index  of  {£„ } is  then  defined  by  the  relation 
H(t)=Hm(t)  (see  Eq.  (1)),  or  more  explicitly 

0(t)=logf/(t)/log//(t),  t E7\ 

Observe  that  0(t)  is  well  defined  since  H has  Uni- 
form[0,l]  marginals  and  hence,  0<//(t)<l  on  T.  The 
following  results  describe  some  basic  properties  of  the 
multivariate  extremal  index. 

Proposition  3. 1 . Assume  that  {f,  } satisfies  A(v„ (t))  for 
each  t GT  and  has  extremal  index  0(t).  Then 

(i)  6(t)=d(tc)  for  each  t GT  and  c>  0,  and 

(ii)  for  each  j=\,...,d,  {£„;}  has  extremal  index  0j=0( t) 
where  t GT  has  all  coordinates  equal  to  I except  the 
jth. 

(Note  that  by  (/),  0(t)  is  constant  along  the  contours 
L,={tc:c>0},  t£7\  and  hence  6j  in  (ii)  is  well-defined.) 

Proof:  Recall  that  (by  Theorem  1.1)  H(tc)=Hft)  and 
H (tc)=Hc(t)  so  that  (/)  follows  from  the  definition  of  the 
extremal  index.  Next,  for  tGT  with  all  coordinates  but 
the  Jth  equal  to  1,  P {Mn^vn(t)}=P  {Mn)<vn](tj)}  and 
hence 

lim  P {Mnj<vnj(tj)}=\im  P {Mn<vn(t)}=H (t)=Hj(tj). 

/i— >oo  n— 

Therefore  by  Theorem  2.2  of  Ref.  [1],  {£„;}  has  ex- 
tremal index  0,  (say)  so  that  Hftf^tf.  Now 
H (t)=Hm(t)  by  definition  of  the  extremal  index,  and  for 
the  present  choice  of  t this  is  the  same  as  Hj(tj)=tfl) , 
whence  it  follows  that  0(t )=9j  for  all  such  t.  □ 

For  tE7\  let  Nf  denote  the  one-dimensional  point 
process  obtained  from  via  the  map  y— from 
{0,1  }d  to  {0,1},  i.e.,  BG®. 

Thus  Nf  has  unit  mass  at  Hn  if  and  only  if  ^,^v„(t). 
Assume  that  {£„ } satisfies  A(v„(t))  and  with  JnX  as  in 
Sec.  2,  let 

iff(y  )=P  {NfUn^y  I Nf(Jn])>0},  y>  1. 

Proposition  3.2.  Assume  that  {£„  } satisfies  A(vn(t ))/or 
each  tG T and  has  extremal  index  0(t).  Then 
0(t)=(linv^2j,ai  yTf^Cy))'1. 


547 


Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


Proof:  Observe  that 


Xv<)0')=£(^,1(-/ni)  I A^‘U,)>0) 


rnP{^v„(  t)> 
(/„.)>  0} 


knrn  nP{^vn(t)} 
n knP{Mrnr£vn( t)> ' 

Now  P {Mn<v„(t)}=//(t)  and 

P {A/„<v„(t)}=//(t)  (by  assumption),  so  that  lim„^ 
«P{£,:£vn(t)}=-log  H( t)  and  (by  Eq.  (13))  linw, 
knP  {Mrn£vn( t)}=-log  H( t).  Therefore  lim,,^ 
yjr^’fy )=l°g  //(t)/log  //(t)=l/0(t),  as  required.  □ 


3.)  While  the  above  results  illustrate  some  of  the  basic 
properties  of  the  multivariate  extremal  index,  they 
are  far  from  complete.  For  instance;  it  would  be 
useful  to  identify  the  set  of  all  “admissible”  O(-)  for 
a given  H , that  is  the  set  of  all  O(-)  such  that  DH(-) 
defined  by  Eq.  (14)  is  a probability  distribution  on 
[0,1  ]''.  It  would  also  be  of  interest  to  study  the  prop- 
erties of  O(-)  when  one  or  both  of  H and  H are 
independent.  In  this  context  we  have  the  following 
simple  result. 

Proposition  3.4.  If  H is  independent,  then  H is  inde- 
pendent if  and  only  if 


Remark:  Proposition  3.2  is  simply  the  multivariate  ver- 
sion of  Eq.  (1)  and  shows  how  the  extremal  index  is 
related  to  the  clustering  of  “exceedances.”  Indeed,  ac- 
cording to  the  present  viewpoint,  an  exceedance  occurs 
at  time  / if  £,^v„(t),  i.e.,  if  ^nj>vnJ(tj)  for  at  least  one  j. 
Thus  Propositions  3.1  and  3.2  show  that  while  the  de- 
gree of  clustering  may  depend  on  t,  it  is  constant  on  each 
Lt.  Note  also  the  connection  to  Theorem  2.8. 

The  next  result  gives  the  relation  between  the  depen- 
dence functions  of  H and  H , which  is  seen  to  involve  the 
extrema!  index  in  an  intricate  manner. 

Proposition  3.3.  //  {£„}  has  extremal  index  0(t),  t £P 
then 

DH(t?\...,tf)=D*i)( t),  ter,  (14) 

where  0j  is  the  extremal  index  of  {£„, },  j=  1 d. 

Proof:  By  definition  of  the  dependence  function, 
T>w(t)-P{//|(Xl)<rl,...y/l/(Xrf)<rrf>  where  (Xu-Xd)  is 
a random  vector  with  distribution  H . Therefore,  since 
HftfMp, 

DH(tfo...,tp)=P{X^t t) 

=DP(t), 

or  required.  □ 

Remarks 

1 . )  Note  that  s=tc  (for  some  c> 0)  if  and  only  if  log  s,/log 

5rf=log  t j /log  f,;=a;  (say),  j=\,...,d—\.  Therefore  we 
may  write  L,=Z.a  where  a=(a,,. and  hence  by 
the  remark  following  Proposition  3.2,  0(t)=0(a), 
i.e.,  the  extremal  index  is  a function  of  d—  1 vari- 
ables only. 

2. )  By  Proposition  3.3,  DH(t^  ,...,tf)=DPt)  (t)=Dfi(tm)- 

Also,  if  t£La  then  (tfo...,tf)GLa,  where  a*=(a,0,l 
0lt,...,aa-\Q,i-\l dj).  Thus  DH  is  obtained  by  translat- 
ing the  values  of  Dh(=H)  on  La  onto  La.. 


d d 

0(t)=^0j log/yX  l°g(/,  for  some  tE.(0,l)d. 

>1 

In  particular,  if  both  H and  H are  independent  then  0(  t) 
is  a convex  combination  of  the  0/s. 

Proof:  If  H is  independent,  then  H (t)=H  (t)m= 
(UPffo'K  The  conclusion  follows  immediately  from 
Corollary  2.7  (iii)  upon  taking  logarithms  and  noting 
that  if  H is  independent,  then  H( t)=II‘l1/l/J. 

The  extremal  index  can  be  given  the  following  more 
general  formulation.  Let  /x  and  p be  the  probability 
measures  on  (0,l)d  corresponding  to  H and  H,  respec- 
tively. Thus  for  instance, 

/x(A)=lim  P{M„Ev„(A)} 

n — 

where  v„(A  )={v„(s)  : s£A  },  AC(0,1)''.  We  now  define 
0(A)  via  the  relationship  p(A)=pe(A)(A),  or  more  di- 
rectly 0(A  )=log/x(A  )/log/2(A ),  for  subsets  AC(0,1)</ 
such  that  fi(A)>0  and  p(A)> 0. 

Note  that  0(t)=0((O,r,)X  — X(0,td))  for  t£P.  Thus  if 
{0(t) : t£P)  is  known  along  with  either  of  H or  H , then 
it  is  possible  at  least  in  theory  to  obtain  {0(A)  : 
ACIO,!/}.  In  practice,  however,  it  may  not  be  possible 
to  obtain  0(A ) in  a tractable  form,  but  frequently  one  is 
only  interested  in  certain  special  sets,  typically  rectan- 
gles of  the  form  Yl^faj  ,bj),  and  for  such  sets  the  compu- 
tation is  easy. 

The  definition  of  Mn  as  the  vector  of  componentwise 
maxima  actually  corresponds  to  regarding  £ as  an  ex- 
treme observation  if  fp-vnj (tj ) for  some  j.  More  gener- 
ally, one  may  define  f to  be  an  extreme  value  if 
fE.v„(A)  for  some  AC(0,l)d,  in  which  case  0(A ) has  an 
interpretation  as  a measure  of  the  clustering  of  such 
extremes.  Note  that  the  original  definition  of  extremes 
corresponds  to  letting  A=((0,r,)X---X(0,rd))c.  Alter- 
nately one  may  consider  taking  A=(q,l)X---X(rrf,l) 
which  corresponds  to  defining  f as  an  extreme  observa- 
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tion  if  €ij>vnj(tj)  for  all  j.  Yet  another  choice  is  A= 
{t : %tf>c  }. 


4.  Examples 


We  conclude  with  two  examples,  both  involving 
bivariate  stationary  sequences. 

Example  4.1  Let  {rj„}  be  an  iid  sequence,  and  put 
£n\=Vn,  and  £„2=max{i7n-i,  }.  Let  F denote  the  distri- 

bution of  £„=(£„\,  £ „i ) with  marginals  F i and  F2.  Then 
F2{x)=P  {£nl<x}=P  {r)n-\<x , rjn<x  }=F2(x)  and 


F(x , , x2)=P  {$,i  <x,  ,£n2—x2  }= 


{ 


F?(x2), 

F{(xi)Fi(x2), 


if  Xj  —X2 
if  X|<r2. 


If  vnj(tj)  satisfies  FJn(vnj(tJ))^>tJ,  7=  1,2,  then  limn_>oo 
F\{vn2{t2))=t22  so  that  vn2(t2)=vn](t22).  Moreover 
vn](ti)>ivn2(t2)  if  and  only  if  tt>t22,  and  so 

H„(vn (t ))=P  {M„<  v„ (t ) }— >"  H(t)=\  2 ’ * 

L ' 1 , 11  <l<<2  • 


The  marginals  of  H are  therefore  // i(/li)=/'i  and  H2(t2)= 
t22  so  that  0i=l  and  02=  1/2,  and  the  dependence  function 
of  H is  DH(t)=H(t\,  tl)=tiAt2.  For  the  associated  iid 
sequence  {£„  } on  the  other  hand,  it  is  easily  verified  that 


Hn (v„ (t))=P  <M„<  v„ (t )={^,„  |J' 

from  which  it  follows  that 


1 


2’ 


loK'l 

log / |/2 


TT 


if  t\>t22 
if  t\<t22. 


We  next  consider  a moving  average  sequence  studied 
in  Ref.  [16]. 


Example  4.2.  Let  {Zk=(Zk]£Zk2)'},  ~^<k<^,  be  a se- 
quence of  iid  random  vectors  in  IR2.  We  assume  the 
existence  of  a sequence  of  positive  constants  an  — >°°,  and 
a measure  v on  IR2  which  is  finite  on  sets  of  the  form 
{x  : ||x||>r  },  r> 0 (where  ||  • ||  denotes  the  Euclidean  norm 
in  IR2),  such  that  nP  {a„  ’ZoE-}— c(-).  (Here  >v’ 
denotes  vague  convergence- of  measures  on  IR2  with 
respect  to  the  metric  d(x\,x2)=  | rj"1  — r2'  | V I 0i~02 1 , 
where  for  i=  1 ,2,r,  and  0,  denote  the  polar  coordinates  of 
x,,  and  a\/b=max{a  ,b  }.)  The  measure  v is  necessarily 
of  the  form  r-({x  : ||x||>r,0(x)E/4  })=r~aS(A ) for  r> 0 
and  AC[0,27t),  where  S(-)  is  a probability  measure  on 
[0,2 7r)  and  a>0.  Hence  in  particular  [17], 


v(cA)=c  av(A)  (15) 

for  all  c>0  and  all  sets  A with  v(A  )<oo. 

Define  the  bivariate  moving  average  process 
X„=X”oC;Z„-7,  where  {C;=[c/W]2>l}^0  is  a sequence  of 
real  2X2  matrices  satisfying  X*/  | Cj#  1 k,l=  1,2,  for 
some  SE(0,a),  §<1.  For  x=(x,,jc2)'EIR2,  write  Axj= 
{z : C,z£((  — oc^xj))1},  where  A1  denotes  the 

complement  of  a set  ACIR2.  Then  [16], 

lim  P {a~'M„<x}=exp{  — y(x)},  and 

/j— »o° 

lim  P {a“'A/n<x}=exp{—  y(x)},  xGIR2, 

n — >°° 

where  y(x)=X”0 v(AXJ)  and  y(x)=^(U”()AXJ).  The  ex- 
tremal index  is  therefore  0(x)=y(x)/y(x),  xEIR2.  It  fol- 
lows from  the  definition  of  AXJ  and  Eq.  (15)  that  this  is 
in  fact  a function  of  xjx2.  Note  that  the  extremal  index 
defined  above  differs  from  that  in  Sec.  3 in  that  it  is 
defined  on  IR2  rather  than  [0,1]2.  However  the  two  defi- 
nitions are  equivalent  as  may  be  seen  by  means  of  a 
suitable  transformation  from  IR2  to  [0,1]2. 

The  actual  calculation  of  0(x)  may  be  quite  difficult 
in  general,  but  possible  to  carry  out  under  appropriate 
simplifying  assumptions. 

Case  (i).  If  C/=CjC  where  C=[cw]2;_,  and  the  c/s  are 
non-negative  constants,  then  AXJ=cj'B(x)  and 
U/,4xj=cB(x),  where  fi(x)={z  : CzE((— °°^Ci)X 
( — °°^r2))c}  and  c=max{cy  : j^O}.  Therefore  by  Eq.  (15), 
v(Axj)=c“v(B(x))  and  v(U^Axj)=cav(B(x))  so  that 
0(x)=c“/X/oc/. 

Case  (ii).  If  the  C/s  are  diagonal,  i.e.,  CJ=diag[c/i,c/2] 
with  Cj,>0,  i=l,2,  then  Axj={z  : cjxZ\>X\  or  cj2z2>x2} 
and  U/<K4XJ=((-°OA'i/ci)X(-co7f2/c2))c  where  c,= 
max{c;,  : 0},  /=  1 ,2.  In  particular,  taking  x2=°°  and 

using  Eq.  (15)  as  in  Case  (i),  we  have  v(AXJ)= 
cj! L’( {z  : Zi>x  1 })  and  v(U^AXJ)=c“v({z  : Zi>x,}),  so 
that  the  extremal  index  of  {X„,}  is  0|=c,7X/oC“,. 

Case  (Hi).  Let  D denote  the  support  of  v.  If  DC 
{z  : z,=0  or  ^2=0}  (which  is  the  case  if  the  coordinates 
of  Z0  are  independent),  then  we  may  write 

c(( - oc^i ) x ( - coyx2)y=a lxl""+<a2x2 ",  x,,x2>0,  (16) 

for  suitable  constants  and  a2> 0.  Once  again,  as- 
suming the  Cjjt/’s  to  be  non-negative  and  writing 
cw=max {cjk!  : y>()}  for  k, 1=1,2,  we  have  (writing 
a/\b=m\n{a  ,b  }) 


AXJC\D=((—,x>rxJ  Cj,\  \ax2!  Cj2\)'Z(~cc^J  Cj_\2/\x2!  Cj_22))cC\D 
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and 

U”0AxjnD=((-'»^i/cii/vc2/c2i) 

X ( - oo^r,/c12Ajr  2/c22))cnD , 

so  that  using  Eq.  (16) 

ai(.x-i/ciiA.v2/c2|)~<'+a2(.Y|/ci2/vv2/c22)-t' 

X ai%*L0(xJCjA\AjC2/Cj'2\)-a+a2Z]Lo(x\/cJ'n*X2/Cj,22ya ' 

Thus  putting  x2= «=,  we  have  6t=(a[c?l+a2c?2)/(alZ'f0cfu+ 
a22“oCjai2)  and  similarly,  02=(o  i c",  +a2c21)l {ax t*LoCf2 , + 

^j-oCf22). 

If  also  Cy,i2=O,2i=0  for  each  j,  (that  is  if  the  C/s  are 
diagonal),  then 

x Ml l£n±£2£T££ 

(3|JCr"S/0c/i  i+<32X2“"X/oC/22  ’ 

and  in  particular,  0i=c“1/X°V‘“n  and  02=c22lXfoCf22.  Note 
that  in  this  case  the  limiting  distributions  of  M„  and  M„ 
are  both  independent,  and  hence  (in  accordance  with 
Proposition  3.4)  9(x)  is  a convex  combination  of  0,  and 
02- 

The  non-negativeness  of  the  C/s  assumed  above  is 
not  crucial  and  may  be  relaxed,  although  at  the  cost  of 
more  involved  calculations. 
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